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ABSTRACT 

In a volume-limited sample of 63 ultracool dwarfs of spectral type M7-M9.5, we have obtained high- 
resolution spectroscopy with UVES at the Very Large Telescope and HIRES at Keck Observatory. 
In this first paper we introduce our volume-complete sample from DENIS and 2MASS targets, and 
we derive radial velocities and space motion. Kinematics of our sample are consistent with the stars 
being predominantly members of the young disk. The kinematic age of the sample is 3.1 Gyr. We find 
that six of our targets show strong Li lines implying that they are brown dwarfs younger than several 
hundred million years. Five of the young brown dwarfs were unrecognized before. Comparing the 
fraction of Li detections to later spectral types, we see a hint of an unexpected local maximum of this 
fraction at spectral type M9. It is not yet clear whether this maximum is due to insufficient statistics, 
or to a combination of physical effects including spectral appearance of young brown dwarfs, Li line 
formation, and the star formation rate at low masses. 

Subject headings: stars: low-mass, brown dwarfs - stars: luminosity function, mass function - stars: 
kinematics 



1. INTRODUCTION 

Ultracool dwarfs - main-sequence stars and brown 
dwarfs of spectral class M7-M9.5 - are of particular im- 
portance for our understanding of Galactic physics for 
a number of reasons. First, M dwarfs are the most nu- 
merous stellar species and represent an important but 
poorly understood contribution to our Galaxy. Our clos- 
est neighorhood is populated by many M dwarfs of which 
a large fraction remains undetected. Second, at least 
three physically important regimes exist across the M 
dwarf spectral classes, just at or close to the tempera- 
ture range of ultracool dwarfs. These are 1) the star- 
brown dwarf boundary; 2) the threshold between par- 
tially and fully convective stars with potential implica- 
tions on the magnetic dynamo and activity; and 3) the 
transition from ionized to predominantly neutral stellar 
atmospheres. We discuss each of these in turn below. 

Stars and brown dwarfs differ by the force that bal- 
ances their gravitational pressure. In a star, central tem- 
peratures and pressure are high enough to start sufficient 
hydrogen burning that stabilizes the star. In a brown 
dwarf, electron degeneracy sets in before hydrogen burn- 
ing so that a brown dwarf cannot reach a stable state but 
becomes cooler and cooler as it ages. Objects less massive 
than M « 0.0 7 are brown dwarfs, m ore massive ob- 
jects are stars (Chabricr & Baraffc 2003). The most mas- 
sive and youngest brown dwarfs can have spectral types 
as early as mid-M but enter the L dwarf regime af ter 
about lGyr ([Burrows et al.lll997t iBaraffe et al.lll998l) . 

In main-sequence stars, the threshold between partially 
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and fully convect ive stars occurs at about M ks 0.3 M© 
(jSiess et al.l l2000h . Partially convective stars harbor a 
tachocline, the boundary layer between the outer con- 
vective envelope and the inner radiative core. The 
tachocline is believed to be the place where the large- 
scale solar, and b y analogy also the sun- like stellar dy- 
namo is situated (jOssendrij verl [2003). This type of dy- 
namo cannot operate in fully convective stars. How- 
ever, no change in magnetic activity or the magnetic 
field strengths is observed among M dwarfs (see, e.g., 
iMohantv fc Basil [2001 iReiners fc Basril 120071 ). On the 
other hand, the timescale for rotational braking and 
the topolo gy of magnetic fields seems to change at this 
threshold ([Reiners fe Basrill2008L [2009] ). 

The third important threshold in the M dwarf regime, 
the transition from ionized to partially neutral atmo- 
spheres, occurs among the ultracool dwarfs of spectral 
type about M9, i.e., at temperatures of about 2400 K 
dMever fe Mever-Hofmeisterl 119991: iFleming et ail 120001 : 
IMohantv et al.ll2002D . There is strong observational ev- 
idence that around spectral type M9, normalized Ha 
emission graduall y decreases with spectral type, i.e. with 
temperature (e.g. IMohantv fe Basril 120031 ) . and that ro- 
tational braking becom es very inefficient at this point 
(Reiners fe Basril[200l . 

This is the first of two papers in which we report on 
high resolution spectroscopy in a volume-complete sam- 
ple of ultracool M dwarfs. In this one, we introduce 
the sample and search for young brown dwarfs via the 
lithium test. We derive radial velocities and report space 
velocities, from which we derive the kinematic age. In a 
second paper, we will investigate the properties related to 
activity, which are rotation, Ha emission, and magnetic 
flux. 
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2. SAMPLE SELECTION AND OBSERVATIONS 
2.1. The sample 

We constructed a sample of known M dwarfs of spec- 
tral type M7-M9.5 that is almost volume-complete (d < 
20 pc). The targets are taken from several catalogues and 
d iscoveries fr o m the DENIS and the 2MASS surveys. 

iCruz et all (|2007l ) present a volume-limited sample 
from 2MASS that contains objects in the spectral range 
M7-L8. Distances are derived from spectro-photometry 
with Mj e stimated from the spectral- type/ Mj calibra- 
tion from ICruz et alJ ([2003). This Mj is combined 
with photometry from the 2MASS Second Incremental 
Data Release PSC to yield Mk s and spectrophotometric 
distances. Uncertainties in distance are dominated b y 
the uncertainties in spectral type (see ICruz 
This sampl e covers ab o ut 40 % of the sky. The sample 
taken from ICruz et al.l (|2007l ) was carefully constructed 
to yield a sample that is complete and volume-limited 
with d < 20 pc in the spect ral range we consider. Nev- 
ertheless, ICruz et all ([2007) note that the sample may 
miss several objects at M7 due to the (J — K s ) > 1.0 
s election crite ri on. 

ICrifo et all (1200! and IPhan-Bao & Bessell (l200fih 
present two parts of a sample from DENIS that covers 
most of the late-M dwarfs up t o spectral type M8.5. Two 
later objects wer e disco vered in lDelfosse et al.l (|2001l) and 
IPhan-Bao et al.l (|2006f ). which make this sample almost 
complete up to M9.5 while covering ~ 13% of the sky. 

For 15 targets, we found parallax measurements in 
the literature, which we use instead of the spectrophoto- 
metric distances where available. One target, 2MASS 
J0041353-562112, was recently reported to show evi- 
dence of a ccretion, whic h implies that it is very young 
(~ 10 Myr; Reiners 2009). At this age, the object is much 
brighter affecting the spectroph otometric dista nce, and 
we use the distance reported by IReinera ([20091 ). Taking 
into account parallax and age information, three of the 
63 targets are found to be located at a distance d > 20 pc. 
We did not exclude these targets from the following in- 
vestigation. 

We constructed a joint sample of objects within 20 pc 
from the two surveys. Both surveys together cover more 
than 50% of the sky and can be considered as almost 
complete in the spectral range M7-M9.5. The full "20pc- 
2MASS-DENIS" survey contains 63 objects, 4 of which 
are found in both the 2MASS and the DENIS surveys. 
The constructed survey probably misses a few M7 and a 
few very late-M dwarfs, and at least three of the targets 
are likely to be located farther away than 20 pc. Never- 
theless, the constructed sample probably provides a rep- 
resentative and robust picture of late-M dwarfs in the 
solar neighborhood. We show the distribution of spec- 
tral types in the sample in Fig.[TJ The sample targets are 
given with their apparent J-band magnitude and spec- 
trophotometric distance in Table[T] 

2.2. Known Binaries 

Our sample contains four binaries that we were previ- 
ously aware of; L HS 1070B (2MASS J0024442-270825B) 
has a companion (|Leinert et al.lll994[) that is about 0.2 " 
away at the time of observation (jSeifahrt et aTl [2008) . 
LP 349-25A ( 2MASS J0027559+22 1932A) was found to 
be a binary bv lForveille et al.1 (|2005[ ). and companions to 



2MASS J1124048+380 805 and 2MASS J2206227-204706 
were discovered by iClose et ah! ([2003D . In all four sys- 
tems, the spectra are dominated by the brighter compo- 
nent so that we carried out our analysis in the same way 
as for single stars. 

During our analysis, we discovered that one of our 
targets, LP 775-31 (2MASS J0435151-160657), shows 
a double-peaked cross-correlation profile indicative of a 
double- lined spectroscopic binary (SB2), i.e., it is a bi- 
nary consisting of two nearly equally bright components. 

2.3. Observations 

Observations for our 63 sample targets were collected 
using the HIRES spectrograph at Keck observatory for 
targets in the northern hemisphere, and using UVES 
at the Very Large Telescope at Paranal observatory for 
targets in the southern hemisphere (PIDs 080.D-0140 
and 081.D-0190). HIRES data were taken with a 1.15" 
slit providing a resolving power of R = 31,000. The 
three HIRES CCDs cover the spectral range from 570 to 
1000 nm in one exposure. UVES observations were car- 
ried out using a setup centered at 830 nm covering the 
spectral range 640-1020 nm. All data provide the Ha line 
as well as the molecular absorption bands of molecular 
FeH around 1 /im that are particularly useful for the de- 
termination of rotation and magnetic fields in M dwarfs. 
These will be investigated in a second paper. 

Data reduction followed standard procedures includ- 
ing bias subtraction, 2D flat-fielding, and wavelength 
calibration using ThAr frames. HIRES data was re- 
duced using the MIDAS echelle environment. For the 
UVES frames, we used the ESO CPL pipeline, version 
4.3.0. There are two important differences to the stan- 
dard pipeline products. First, the version of the pipeline 
used in Period 81 produced wavelengths that are offset 
from the correct solution by a constant factor of 1 A, this 
bug has been fixed in version 4.3.0. Second, the standard 
pipeline extracts the spectra before flat-fielding, which is 
not appropriate for spectra with strong fringing patterns. 
The new pipeline version allows to flat-field the original 
2D frames so that fringe patterns are effectively removed. 

The apparent brightness of our targets arc in the range 
J = 8.9-13.4mag. Exposure times between 200s and 3h 
yielded signal-to-noise ratios between 20 and 100 at 1 fim 
in all objects. 

3. MEASUREMENTS 
3.1. Li Absorption 

Our spectra cover the 6708 A Li line that can be 
used for the "lithium-test" in l ow-mass objects (e.g., 
IMaeazzu et al.l 119911: IBasril I2000h . In our spectra, the 
detection limit of the Li line is around 0.5 A. While stars 
quickly deplete Li on a timescale of 100 Myr and less, 
brown dwarfs need longer for this process or do not en- 
tirely _depJete_IJ_a^_alL__We show evolutionary tracks 
from lD'Antona fc Mazzitellil (|1997D in Fig.H The dashed 
blue line marks the region where Li is depleted by 99 %, 
i.e., objects with detected Li lie on the left hand side 
of this line. For the temperatures of our sample stars, 
2400 K< T off < 2800 K, this means that objects with Li 
are less massive than about 65 Mj , and younger than 
~ 0.5 Gyr. This also means that they are less massive 
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than the substellar limit at 75 Mj , hence they are young 
brown dwarfs. 

We have detected Li in six of the 63 objects in our 
sample, i.e., about 10% of our sample are brown dwarfs 
younger than about 0.5 Gyr. The Li lines of all six 
objects are shown in Fig.[3J Only one of the targets, 
LP 944-20 (2MASS J0339352-352544), was previously 
known as a y oung brown d warf through the detection 
of its Li line (|Tinnevl II 998). Youth and brown dwarf 
nature of the oth er five objects were unknown before. 
ICruz et all (|200l classified 2MASS J0443376+000205 as 
a low gravity, hence probably young, object. This classi- 
fication is supported by our Li detection. 

3.2. Radial Velocities and Space Motion 

Radial velocitites were measured relative to a spectrum 
of a standard star. We used Gl 406 as radial veloc- 
ity standard (w rad = 19 ± lkms -1 , IMartm et all 119971) 
using the cross-correlation technique and correcting for 
barycentric motion. We corrected for wavelength calibra- 
tion inconsistencies by cross-correlating telluric lines in 
the O2 A-band. Thus, uncertainties in our radial veloc- 
ities are on the order of 1-3 km s -1 depending on rota- 
tional line broadening. The uncertainties of the spec- 
trophotometric distance are typically on the order of 
~ 10 %. For our sample of nearby objects this introduces 
an uncertainty in the space motion components of typ- 
ically ~ lkms -1 . If an object is substantially younger 
than assumed for the spectrophotometric distance calcu- 
lation, the error is much larger. This effe c t is i llustrated 
for 2MASSJ 0041353-562112 in IReinersI (f2009h . but is 
probably not important in the other objects. 

Space motions U, V, and W were computed using the 
IDL procedure gal_uvw 2 , but we use a right-handed 
coordinate system with U towards the Galactic cen- 
ter. For most of our targets , proper motions are re - 
port ed inlFahertv et all (120091) . |Phan-Bao et alj (|2001l ) 
and iPhan-Bao et alj (|2003l ) measured proper motions 
for some of our objects, and we use their values where 
available. One object, LHS 1070B (2MASS J0024442- 
270825B), is not included in eith er of the works, and w e 
use the proper motion given in Salim fe Gouldl (2003). 
All radial velocities and space motion vectors [U, V, W] 
together with the references for distance and proper mo- 
tion are given in Table[TJ 

4. RESULTS 
4.1. Lithium Brown Dwarfs 

From our volume-limited sample, we can draw robust 
constraints on the fraction of young brown dwarfs among 
the ultra cool dwarfs . Our results are summarized in 
Table[2] IReid et all (|2002af ) found fractions of F hi = 
(6 ± 4)%) for M7-M9.5 dwarfs and F Li = (10 ± 7)% for 
M8-M9.5 dwarfs. The Li fractions from our sample are 
somewh at higher than the ones reported in IReid et aLl 
(|2002aD . but the discrepancy does not exceed the uncer- 
tainties due to the small sample size. 

We conclude that the fraction of lithium brown dwarfs 
among M7-M9.5 dwarfs is at most only little higher than 
previously reported. This means that the mass function 
index a could perhaps be a little higher than reported in 

2 http://idlastro.gsfc.nasa.gov/contents .html 



IReid et all (|2002af ). Nevertheless, their conclusion that 
a < 2 for both the A r izona and Lyon models is still valid. 

iKirkpatrick et all (|2008f) investigated the fraction of 
objects with Li detections among L dwarfs together with 
a search for objects with spectroscopic signatures of 
youth, i.e., low gravity. They show that the fraction of 
Li detections is about 10 % at spectral type L0 and rises 
to about 50-60% at mid-L (see their Fig. 17). After L6, 
the fraction turns over and declines towards the late- 
L types proba b ly because Li disappears into molecules 
(|Lodd crs 1999). IKirkpatrick et all ((2008) emphasize that 
Li is probably not detectable in very young L dwarfs, i.e., 
during the first few ten million years, because the Li line 
is not visible at very low gravity. 

We plot the fraction of lithium brown dwarfs as a func- 
tion of spectral types M7- L8 in Fig.S ] M7-M 9 dwarfs 
are from this w ork and from | Reid et al.l ([2002a) , L dwarfs 
are taken fromlReiners fe Basril (l2008h . and from the Li 
analysis of IKirkpatrick et alj (I2008D including data from 
IKirkpatrick et all (|1999ll2000D . While the Kirkpatrick et 
al. works employ low-resolution LRIS spectra, the other 
Li detections are based on spectra with higher spectral 
resolution. Nevertheless, the detection thresholds from 
the different samples are not very different because in L 
dwarfs it is usually not the spectral resolution that limits 
the detection of Li (note that the Li lines shown in Fig. [3] 
are smoothed to a lower effective resolution). Compar- 
ing the subsamples taken from high- and low-resolution 
spectra, we see no difference in the fraction of Li detec- 
tions, and the results for objects that are contained in 
both samples are consistent. 

The fraction of lithium brown dwarfs among ultracool 
dwarfs continually rises from M7 to M9, because late-M 
dwarfs in general are less massive, so that the fraction of 
objects with M < 0.07 M Q is growing towards later spec- 
tral types. In Fig.[H there is also a general trend that the 
Li fraction rises from M7 to mid-L, which is consistent 
with a higher brown dwarf fraction at later spectral type. 
Note that the time required to deplete Li in high-mass 
brown dwarfs (> 60 Mj) is longer among the L dwarfs 
(~ lGyr) than at late-M spectral type (~ 300 Myr). 
However, the fraction of Li detections among M9 dwarfs 
is much higher than at M8, and, in particular, it is much 
higher than at spectral types L0 and LI. The number of 
objects in this sample is still fairly low, and the results 
are formally still consistent with a smooth rise in the 
range M7-L6 if the M9 bin is interpreted as an outlier. 
On the other hand, the discontinuity at spectral type M9 
is fairly high while the overall scatter between the other 
bins is rather small. Three out of the 15 objects at spec- 
tral type M9 and M9.5 show Li, all three are most likely 
within 20 pc. The statistical significance of the M9 dicon- 
tinuity is comparable to that of the apparent downtur n 
in Li at late L types noted by IKirkpatrick et all ((2008) . 
In the following, we discuss possible explanations. 

The first possible explanation is that our combined 
sample is not strictly volume-limited, in particular the 
L dwarfs are not selected according to their distance. 
This implies a bias towards young objects and proba- 
bly causes an excess of young dwarfs that still have Li. 
Most distances of our M dwarf sample are taken from 
spectrophotometry, which means that a few other young 
objects contained in our M dwarf sample may be situ- 
ated outside 20 pc. However, this problem should apply 
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even more to the early-L dwarfs, in which no effort was 
made to include faint old objects out to a certain dis- 
tance. Thus, one would expect the early-L dwarfs to 
be biased towards higher fraction of Li detections, and 
this bias should be higher than the bias among the late- 
M dwarfs. Furthermore, the fraction of Li detections in 
our M9/M9.5 targets from the volume- limited sample is 
higher than the fraction in the full sample including the 
targets from lReid et al.l (|2002af ). Thus, an observational 
effect that causes a bias towards young objects particu- 
larly at the M9 spectral bin is unlikely an explanation 
for the observed discontinuity. 

A second potential reason for the discontinuity is that 
the L dwarf spectra are of lower quality compared to 
the M dwarf spectra taken for this work. This would 
lead to a higher Li detection threshold and hence a 
lower detection rate. While such an explanation might 
hold for the low-resolution spectra (a detection thresh- 
oldof 3 A is reported in iKirkpatrick et al.l 119991 120001 . 
I2008D. it does not apply to the high resolution spectra 
of I Reiners &: Basrl (|2008f ) . In the latter sample alone, 
the ratio of Li detections to total number of targets is 
0/10 (0%) at spectral type L0 and 1/14 (7%) at spec- 
tral type LI, which is an even lower ratio than t he on e 
from low resolution spectra in IKirkpatrick et all (2008). 
Thus, differences in data quality do probably not lead 
to systematic differences in the Li fraction, and it seems 
unlikely that the detection threshold is the main reason 
for the low fraction of Li detections at early L dwarfs. 

As a third potential reason, we speculate that the dis- 
continuity is real, i.e., the fraction of objects with de- 
tectable Li peaks at spectral type M9. We see two possi- 
ble ways to explain this. 1) Lithium brown dwarfs with 
temperatures of ~2200K may be classified as M9, but at 
older ages stars of the same temperature could be clas- 
sified as L0. Classification of spectral types between M9 
and L0 is based on molecular species that are temper- 
ature dependent and sensitive to the presence of dust, 
which is strongly influenced by gravity. In other words, 
the relation between temperature and spectral class may 
be a function of age. 2) The history of the star formation 
rate toget her with the detectabili ty of Li as a function 
of age (see IKirkpatrick et al.ll2008f ) may conspire so that 
today we observe an excess of Li detections at spectral 
type M9. This may imply a burst of star formation at 
a particular age. So far, however, we are not convinced 
that either of the two is a very likely scenario, and we 
will have to wait for a larger sample and better distance 
estimates to see whether the maximum disappears with 
better statistics. 

4.2. Space Motion of Ultracool Dwarfs 

Space motions of M dwarfs and ultracool dwarfs have 
raised a lot of interest during the last decades because 
the velocity dispersion of a sample can be tied to its age, 
which allows the study of the star formation rate and the 
time-dep endence of physical processes like for example 
activity. lHawley et al.l (|1996l ) have shown in a volume- 
limited sample of M dwarfs that the velocity dispersion 
of active dMe sta rs is smaller than th e dispersion of older 
non-active stars. iReid et all (|2002af) determined the ve- 
locity dispersion of ultracool objects finding that the ve- 
locity dispersion of their M7-M9.5 sample is smaller than 



the d isp ersion in t h e full dM sample of lHawlev et al.l 
(|1996f) . IReid et al.l (|2002bD revisit the space velocities 
of ultracool M dwarfs suggesting that the early-M dwarf 
sample consists of members of both the thick disk and 
the thin disk. Late-M dwarfs, on the other hand, may 
have lost their thick disk component because the lower 
metallicity of these stars moves the hydrogen burning 
limit to earlier spectral types so that these stars have 
cooled towards later spectral types and don't show up as 
M7-M9.5 objects any longer. 

We show the space velocities in U — V and U — W 
diagrams in Fig.[5j and the distribution of U, V, and W 
in Fig. [51 Lithium brown dwarfs are shown as red stars 
in Fig. [5] and are overplotted as filled red histograms in 
Fig. El The distribution of space velocities in ultracool 
dwarfs is clearly concentrated around the thin disk, and 
the six lithium brown dwarfs cluster in a very narrow 
range in the center of the thin disk in all three space 
velocities. This clustering implies a very small velocity 
dispersion, w hich is an indicator for a young sample (see 
iWielenl [19771 and next Section). Thus, space velocities 
of lithium brown dwarfs are consistent with very low age. 
In the following, we will derive the age of our sample of 
M7-M9.5 dwarfs. 

4.3. Ages from Velocity Dispersions 

IWielenl (|1977l ) has shown how a velocity dispersion of 
a sample of stars can be used to determine t he age of the 
sample. More recently, iFuchs et al.l (|2001h refined this 
picture using more recent data from different works. 

First, i t is important to realize that the formulae pro- 
vided by Wielenl (|1977f) are valid for dispersions ajj, cry, 
and aw that are weighted by their vertical velocity \W\ 
(Eqs. 1-3 in IWielenl 1 19771 see also Reid et al., 2002b). 
Although it might be arguable whether this weighting 
is useful or not when comparing dispersions from differ- 
ent samples, it has to be applied if the age of a sam- 
ple sh ould be determined from the formulae in IWielenl 
(1977). The disper sion of space velocities in our sam- 
ple is (ay, oy, aw) — (27.2, 22.6, 14.6) krns" 1 , and the 
| W | -corrected velocity dispersions are (au, cry , o~w) — 
(30.7, 22.6, 16.0)kms" 1 . 

A second but much more crucial principle for the calcu- 
lation of ages from velocity dispersions is the calculation 
of the total velocity dispersion, 

<7tot = {ou + °v+Ow) 1/a - (!) 

This means that <7 to t is not the dispersion of the total 
velocity scalar u tot = iJJ 2 + V 2 + W 2 ) 1 ^ 2 , since the lat- 
ter will always be smaller because motions into different 
directions partially cancel out each other (for example, a 
sample with velocities equally distributed into all direc- 
tions but with the same total velocity will have zero dis- 
persion around Vtot)- The uncorrected total dispersion of 
our sample is crtot = 38.3kms _1 , and the |VF|-corrected 
total dispersion is otot = 41.3 km s . 

It should be noted that, following the same argumen- 
tation as for <7tot, the total dispersion in the tangential 
velocity scalar, i>t a n, cannot be applied to calculate otot 
using otot = (3/2) 1 / 2 (Ttan: The total dispersion in the 
two-dimensional velocity v tan must be the squared sum 
of the two individual dispersion components. Thus, the 
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total dispersion in u ta n cannot be used for the calculation 
of age from the Wielen-relation. 

Recently, kinematics of ultrac ool dwarfs and brow n 
dwarfs have been pre sented by ISchmidt et al.l ([20071. 
IZapatero-Osorio et al.l (|2007D . and iFahertv et all (|2009h 
using the velocity dispersion of a sample of stars to 
cal culate the age o f the sample fo llowin g the work 
of iWielenl (|1977T ). ISchmidt et all (|2007h derives an 
age of 3.1Gyrs for a s ample of late-M and L dwarfs, 
IZapatero-Osorio et alJ (|2007f ) finds ages of 4Gyrs for 
thei r late-M dwarfs and a bout 1 Gyr for L and T dwarfs, 
and [Faherty et al. (2009) reports ages between 2 and 
3 Gyrs for late-M, L, and T dwarfs. Unfortunately all 
three of these papers contain the errors in interpretation 
of the total velocity dispersion pointed out in the previ- 
ous two paragraphs. 

Once the velocity dispersion otot 1S determined for a 
sample , it ca n be translated into an age according to 
IWielenl (|1977l ). He presents three different equations to 
carry out this task, and it is often difficult to find out 
which prescription authors use when publishing an age 
for a sample of stars. The first prescription assumes a 
constant diffusion coefficient an d is apparentl y not much 
used in recent work (Eq. (8) in IWielenl [19771. The two 
other formulae assume velocity-dependent diffusio n coef- 
ficient s. The two formulae, Eqs. (13) and (16) in IWielenl 
(fl977l) . differ in the way they treat the diffusion coeffi- 
cient C v , they are 

3 

^ = <o+2^ T ( 2 ) 

4 = <o + llvMexp (r/T 7 ) - 1), (3) 

with t the age of the sample in yr, a v .o = lOkms -1 , j v = 
1.4 10" 5 (km s- x ) 3 /yr, >y ViP = 1. 1 10~ 5 (km s" 1 ) 3 ^, 
and T 7 = 5 10 9 yr. According to IWielenl (fl977f ). Eq.[2] 
is inadequate for ages above 3 Gyrs, which means that 
Eq. [3] should be preferred 3 . 

The determination of an age from a stellar sample as- 
sumes that the sample members (at least in a very qual- 
itative sense) fulfill the requirement of homogeneity, i.e., 
their velocity distribution should not differ drastically 
from a Gaussian distribution. A robust way to investi- 
gate the velocity distribution in a sample is to use p roba- 
bility plots, or "probit" plots (jLutz fc Upgrenfl "l980). In a 
probability plot, the data points are sorted and assigned 
a probability bin, which is the distance to the mean in 
units of the standard deviation, that each of the sorted 
points in a strictly Gaussian distribution would have. If 
a distribution is Gaussian, the sample points will follow 
a straight line, and the value of the standard deviation, 
er, is the steepness of this line. We show the probabil- 
ity plots for our sample in Fig. [7] The largest part of 
the sample relatively well resembles a well-constrained 

3 ISchmidt et al.l H2007I) use Eq. (0 and present it in a slightly 
different form with t the age of the sample and r a characteristic 
timescale, r = 2 10 8 yr. However, this actually should have been 
written as r~ 1 rj 2 10 _8 yr _1 . It provides a fortunate compen- 
sation for th e error in comp uting tot al velocity dispersions. Thu s 
the results of | Schmidt et alJ 1120071 '). IZapatero-Osorio et ail (|2007l) . 
and lFahertv et al.l ( 1200917 end up fortuitously c ons istent with other 
work (and ours). To illustrate, Schmidt ct al. (2007) find an age of 
3.1 Gyr for otot = 20.8kms _1 , but this low a velocity dispersion 
would in fact imply an age well below 1 Gyr. 



line, which means that except for a only a handful of 
outliers, the distribution of space velocities is consistent 
with a Gaussian distribution. We overplot a fit to the 
inner ±1.75<r of the velocity distribution in Fig. [71 Ve- 
locity dispersions and ages derived from the inner part of 
this distribution are statistically indistinguishable from 
results taken from the full sample. 

From the \W\ -corrected space velocity dispersions we 
have calculated the age of our sample as described above. 
Our results are summarized in Tabled the age estimate 
for our volume-limited sample of ultracool dwarfs from 
this methodology is r = 3.1 Gyr. We compare thi s 
result to the sample presented in iReid et alT ((2002a) . 
These authors also show that their sample matches a 
Gaussian distribution quite well. We calculated |W|- 
corrected space velocity dispersions from the 32 M7- 
M9.5 dwarfs (cp. their Eq. (3)) given in their Table 5, 
and from the Wielen-relation we derive an age for the 
sample of t = 3.1 Gyr. Given the substantial uncertain- 
ties of this age-determination, which are on the order 
of several hundred Myr, this is in remarkable agreement 
with the age derived from our sample. Thus, there is 
good evidence that the mean age of the local ultracool 
dwarf population is about 3 Gyr. 

5. SUMMARY 

Ultracool dwarfs of spectral types M7-M9.5 are of par- 
ticular interest to our understanding of the physics of 
low-mass stars and brown dwarfs. In order to investigate 
the physics of ultracool dwarfs, we constructed a volume- 
limited sample of 63 M7-M9.5 dwarfs from the DENIS 
and 2MASS samples. In this first paper, we present the 
sample composition, the fraction of young brown dwarfs 
from detections of Li, and the kinematic age of the local 
population of ultracool dwarfs. 

We detected the Li line in six of our targets, i.e., about 
10% of our targets are younger than about 0.5 Gyr. This 
result is consisten t with ear li er resu lts on the brown 
dwarf fraction by IReid et all (|2002afh The mass func- 
tion index a is perhaps a little higher than previously 
thought, but a < 2 is probably still valid. 

From M7 to mid-L, the Li fraction rises because the 
cooler stars are less massive and retain their Li for a 
longer period. However, we find a local maximum of the 
Li fraction at spectral type M9. At this point, we cannot 
exclude that this maximum is caused by the full sample 
being partially biased towards younger objects that are 
brighter and hence easier to observe. A potential expla- 
nation for a real Li maximum could be a situation in 
which the spectral appearance of lithium brown dwarfs 
together with the mass-dependent detectability of the Li 
line and the star formation rate results in a maximum of 
stars with detectable Li lines at spectral type M9. The 
details of such a scenario are not clear, and a statistically 
more robust sample of late-M and early-L dwarfs should 
be investigated before drawing further conclusions. A 
necessary condition for this is the availability of trigono- 
metric parallaxes. 

We investigated the space velocities of ultracool 
dwarfs. UVW- velocities are concentrated around the 
thin disk, and the lithium brown dwarfs kinematically 
cluster around the very center of the thin disk. Space 
velocities of our sample generally follow a Gaussian dis- 
tribution so that we can derive a kinematic age for our 
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sample. The kinematic age of M7 -M9.5 dwarfs is 3.1 Gyr, 
which supports earlier results (see lReid et al.ll2002al) that 
the velocity dispersion of late-M dwarfs is significantly 
lower than those measured for nearby M dwarfs, and 
that the dispersion is similar to the one in earlier type 
emission-line stars. 

Based on observations collected at the European 
Southern Observatory, Paranal, Chile, PIDs 080.D-0140 
and 081.D-0190, and observed from the W.M. Keck Ob- 
servatory, which is operated as a scientific partnership 
among the California Institute of Technology, the Uni- 
versity of California and the National Aeronautics and 
Space Administration. We would like to acknowledge 
the great cultural significance of Mauna Kea for native 
Hawaiians and express our gratitude for permission to 
observe from atop this mountain. We thank Andreas 
Seifahrt and Denis Shulyak for helpful discussions on 
space velocities and Li abundances. A.R. has received 
research funding from the DFG as an Emmy Noether 
fellow (RE 1664/4-1). G.B. thanks the NSF for grant 
support through AST06-06748. 
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Sample targets and their space motion. 
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21 


-64 


-9 


-61.2 




2MASSJ 


2306292 




050227 




M8.0 


11.35 


11.00 


-41 


-64 


19 


-56.3 




2MASSJ 


2331217 




274949 




M8.5 


11.65 


11.60 


-22 


36 


3 


-4.2 




2MASSJ 


2349489 


+ 


122438 


LP 523- 55 


M8.0 


12.60 


19.60 


8 


-15 


-10 


-3.5 




2MASSJ 


2351504 




253736 




M8.0 


12.47 


18.20 


-37 


-3 


3 


-10.0 




2MASSJ 


2353594 




083331 




M8.5 


13.03 


19.10 


49 


36 


-18 


32.7 





Ref eren ces. — fal IReinersI (2QQ3), (b) ICosta et all pOOST ). Cc) ICosta et al.l |[2lX)6T) . (d) ILepine et all (t2009l'). (e) IMonet et al.l 
fl99^l. ({) IBartletti J20Qg) ; proper motion usual ly taken from IFahertv et all ( 120091) , except for: fl) IPhan-Bao et al.l (120011 ). (2"T 
IPhan-Bao et all l |2003f l. f3 nSalim fc Gouldl pOOST ) 

Note. — Li: lithium brown dwarf - SB1: single lined binary - SB2: double lined binary 

* The typical uncertainty for ti ra< j is 1— 3kms — 1 depending on rotational line broadening. 
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Fig. 1. — Spectral type distribution of the sample targets. 
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TABLE 2 
Li fractions 



Spectral bin 


N u 


iVtot 


Li fraction [%] 


M7 - M9.5 


6 


63 


io±t 


M8 - M9.5 


5 


39 


13±l 


M7/M7.5 


1 


24 


4 +9 


M8/M8.5 


2 


27 


7+9 
'-2 


M9/M9.5 


3 


12 


25l« 



3000 



2900 - 



2800 r 



2700 - 



2600 - 



2500 - 



2400 




M6 



M7 



_ M8 



M9 



50 



100 



150 



250 



300 



200 
Age [Myr] 

Fig. 2. — Evolutionary tracks of low mass stars and brown dwarfs from D'Antona & Mazzitclli (1997), we chose the 1998 updated models 
with a Deuterium abundance Xn = 10~ 5 (Ccnsori & D'Antona 1998). Tracks are shown in red where the Li abundance is depleted by 
99% and more. The dashed blue line divides the plot into regions where L i is not yet depleted (left of the dashed line) and it is depleted 
(right). Approximate spectral types according to Golimowski ct al. (2004) are indicated on the right axis. Note that this plot covers the 
mid- and late-M dwarfs and that the L dwarf regime is not shown. 
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6700 6705 6710 6715 
Wavelength [A] 



6700 6705 6710 6715 
Wavelength [A] 



6700 6705 6710 671 
Wavelength [A] 



Fig. 3. — Li lines, from left to right: top row: 2MASS J0019262+461407 (M8.0, note that rapid rotation causes the broad and relatively 
shallow Li line), 2MASS J0041353-562112 (M7.5)), 2MASS J0123112-692138 (M8.0); bottom row: 2MASS J0339352-352544 (M9.0), 
2MASS J0443376+000205 (M9.0), 2MASS J1411213-211950 (M9.0). Spectra have been smoothed with a 15-pixcl box-car for better 
visibility. 
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Fig. 4. — Fraction of objects with detected Li per spectral type. M7-M 9 objects are from this work and Rcid ct al. (2002a), L dwarfs 
are from [Rcincrs & Basri (2008) and Kirkpatrick ct al. (1999, 2000, 2008). Numbers denote the total number of objects per spectral bin, 
error bars are la uncertainties. 
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TABLE 3 

Velocity dispersions au,crv, and aw, total dispersion 

ftoti AND AGE CONSTRAINT S FROM THIS S AMPLE AND FROM 

the M7-M9.5 sample from IReid et al.I (I2002AI ) . Velocit y 

DISPERSIONS ARE |H^|-CQ RRECTE D ACCORDING TO IWlELENi 
lfT977l) . 







°V 


a w 


CTtot a 


Age 




[kms -1 ] 


[kms -1 ] 


[kms -1 ] 


[kms -1 ] 


[Gyr] 


this work 


30.7 


22.6 


16.0 


41.3 


3.1 


Reid et al. (2002a 1 ) 


33.7 


17.4 


17.4 


41.7 


3.1 



a atot =(^+4 + 4)^ 



P V c 
o^_*S-^f ° o 

° «P 

O O o 




100 



50 



-50 




-IOOLl- 
-100 



-50 





[km/s] 



50 



100 



■IOOLl. 
-100 



-50 





[km/s] 



50 



100 



Fig. 5. — Space velocities in U — V (left) and U — W (right) diagrams. The thin disk is indicated as a red ellipsoid, and the six lithium 
brown dwarfs are shown as red stars. 
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Fig. 6. — From left to right: Histograms of space motion components U, V, and W. The distribution of lithium brown dwarfs is overplottcd 
as a filled red histogram. 
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